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With this paper we propose a model to simulate the functional aspects of light adaptation in 
retinal photoreceptors. Our model, however, does not link specific stages to the detailed molecular 
processes which are thought to mediate adaptation in real photoreceptors. We rather model the 
photoreceptor as a self-adjusting integration device, which adds up properly amplified luminance 
signals. The integration process and the amplification obey a switching behavior that acts to 
locally shut down the integration process in dependence on the internal state of the receptor. The 
mathematical structure of our model is quite simple, and its computational complexity is quite low. 
We present results of computer simulations which demonstrate that our model adapts properly to 
at least four orders of input magnitude. 
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I. INTRODUCTION 

There is agreement that adaptation (i.e. the ad- 
justment of sensitivity) is important for the function 
of nervous systems, since without corresponding 
mechanisms any neuron with its limited dynamic 
range would stay silent or operate in saturation 
most of the time [jj. Because neurons are noisy 
devices, reliable information transmission is only 
granted if the distribution of levels in the stimulus 
matches the neuron's reliable operation range 0. 
Consider, for example, the mammalian visual 
system, with the retina at its front-end. When 
performing saccades, the retina must cope with 
intensity variations which may span about one 0, Q 
to about two orders of magnitude (2 including shad- 
ows acoording to 0, 2-3 according to [j|)- From 
one scene to another (e.g., from bright sunlight to 
starlight), the range of intensity variations may well 
span up to ten orders of magnitude |fj, [D, |8, 9] ) . This 
range of intensities has to be mapped onto less than 
two orders of output activity of retinal ganglion 
cells |lOj , implying some form of compression of the 
scale of intensity values. The retina achieves this by 
making use of a cascade of gain control and adapta- 
tion mechanisms, respectively (e.g., [ill Ha . Ha. Ill ] ) . 
Specifically, cone photoreceptors may decrease their 
sensitivity proportionally to background intensity, 
over about 8 log units of background intensity . 
This relationship is known as Weber's law (e.g., 
fll|). Adaptation in photoreceptors is achieved 



by subtly balanced network of molecular processes 
(see 01 f° r an excellent introduction, and [T3. ITsf 
with references). 

With the present paper, we propose a mechanism 
which mimics the dark and light adaptation of 
retinal cones. Our mechanism abstracts from the 
detailed molecular processes of the transduction 
cascade as described in the following section. We 
scckcd out an easy implementable and computa- 
tionally efficient way of achieving the adaptation 
behavior of cone photoreceptors. Our approach 
should - and will be - contrasted with the retinal 
stage of a recently proposed model of Grossberg 
and Hong 0, |2(j, which simulates (i) luminance 
adaptation at the outer segment of the photore- 
ceptor (c.f. |2l|). and (ii) inhibition at the inner 
segment of the photoreceptor by horizontal cells 
(e.g., [12 )■ In their model, horizontal cells are 
coupled with gap junctions (forming a syncytium), 
whose connectivity or permeability decreases with 
increasin g di fferences between the input of adjacent 
cells [23l 124). In other words, their horizontal cell 
network establishes current flows inside of regions 
that are defined by low contrasts, whereas no 
activity exchange occurs between regions which are 
separated by high contrast boundaries (very similar 
to an anisotropic diffusion mechanism |25|). In this 
way, contrast adaption is implemented. Notice that 
our model lacks the latter stage, and only simulates 
the photoreceptor adaptation. 



"This manuscript (co-authored with Jordi Vitria) has been 
submitted to EURASIP Journal on Applied Signal Processing 
- Special Issue on Image Perception, to appear 3rd Quarter 
2006; Electronic address: threequarks@yahoo.com 



II. MECHANISMS OF ADAPTATION IN 
THE RETINA 

A response to light is initiated by photoiso- 
merization of the chromophore ll-cis-retinal to 
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alKrcms-retinal. In darkness, 11-cis-retinal is 
bound to rhodopsin in its inactive conformation, 
and lies buried in the membranes of the outer seg- 
ment discs. Upon absorption of a photon, and the 
subsequent photoisomerization of the chromophore, 
the rhodopsin undergoes a conformational change 
which converts it into its active form Rh* (or 
metarhodopsin II). The presence of Rh* triggers 
two distinct mechanisms: an recycling process 
known as visual cycle, and an enzymatic cascade 
known as transduction cascade. 
The visual cycle begins with the phosphorylation of 
Rh* , and subsequent binding of arrestin to the phos- 
phorylated photopigment. After binding of arrestin, 
the photopigment is rendered completely inactive. 
The protein opsin is then dephosphorylatcd, and 
alKrcms-retinal is reduced to all-trcms-retinol. The 
retinol is isomerized to the 11-cis-isomer outside 
the photoreceptor (in the adjacent retinal pigment 
epithelium layer), and re-enters afterwards to 
re-combine with the dephosphorylated opsin. 
The transduction cascade begins with the serial 
activation of transducins by Rh* , implementing the 
first stage for signal amplification in the cascade 
[2bT |. Thereby, an active complex T a -GTP is 
formed, which binds to and activates the enzyme 
phosphodiesterase (PDE). PDE reduces the con- 
centration of cytoplasmatic cGMP by hydrolizing 
it. The latter process constitutes a second stage 
for amplification. The hydrolysis of cGMP causes 
the closing of cGMP-gated channels, what in turn 
generates the electrical response of photoreceptors. 
Thus, photoreceptors are depolarized in darkness 
because of their open cationic channels, and get 
hyperpolarized by light. In darkness, the 

steady current that flows into the outer segment 
is usually called dark- or circulating current |4lj . 
The main fraction of the circulating current is 
carried by Na + ions, and a smaller fraction of Ca 2+ 
ions [2]j. Calcium is transported out of the outer 
segment by the Na + /K + -Ca 2+ -exchange protein 
at a constant rate, independent of the light hitting 
the photoreceptor. This implies that light decreases 
intracellular Ca 2+ levels, because of the increased 
probability of channel opening. As a consequence, 
a direct correlation (i.e., a linear relationship) 
exists between the circulating current and Ca 2+ 
concentration. 

Adaptation of the photoreceptor to ambient light is 
granted by balancing the just described amplifica- 
tion mechanisms (for low light situations) against 
mechanisms which prevent response saturation (e.g., 
for sunlit scenes). This balance is implemented 
by feedback mechanisms which act either on the 
catalytic activity or on the catalytic lifetime of the 
components that make up the phototransduction 



mechanism 


Ref.[19, 20} 


our approach 


light adaptation 


yes 


yes 


divisive gain control 


yes 


no 



TABLE I: Model overview. The table gives an 
overview over the mechanisms used in the model of G&H 
[l9l l2(i|l and our approach. 



cascade |28| • It is now well established that changes 
in Ca 2+ concentration regulate the cascade in at 
least three important ways. 

First, Ca 2+ can prolong the lifetime of Rh* through 
the inhibition of phosphorylation in the visual cycle, 
by means of recoverin. Second, in the transduction 
cascade, Ca 2+ regulates the cytoplasmatic concen- 
tration of cGMP by binding to guanylate cyclase - 
the enzyme that is responsible for cGMP synthesis. 
Third, decreasing Ca 2+ concentrations increase the 
sensitivity of the cationic channels to cGMP [2jJ . 
Taken together, Ca 2+ is now considered as the 
photoreceptor's internal messenger for adaptation. 
Supporting evidence comes from the fact that 
adaptation effects can be provoked without light 
(c.f. p. 130 in 01 )i by only lowering the Ca 2+ 
concentration, or that adaptation is suspended by 
clamping the Ca 2+ level to its value corresponding 
to darkness (p. 126 in [b^b 

Beyond the level of the individual photoreceptor, 
further mechanisms related to adaptation are 
effective, for example network adaptation in in- 
terneurons and retinal ganglion cells (i.e. adaptation 
is "transferred" beyond the rece ptiv e field of the 
actually stimulated cell, e.g. |30l l3ll I32I l33j "> . and 
discounting predictable spatio-temporal structures 
from the stimulus by Hebbian mechanisms [34l l35| . 



III. FORMAL DEFINITION OF THE 
ADAPTATION DYNAMICS 

Table [I] gives a brief comparison of components, 
and a sketch of our model is shown in figure ^ I n 
what follows, we give the formal introduction to our 
mechanisms which are thought to provide an ab- 
stract view for adaptation as it takes place in the 
outer segment of individual photoreceptors. 
Let Cij be a two-dimensional luminance distribution 
which provides the input into our model. For the 
purpose of the present paper we assume that the 
model converges before changes in luminance occur, 
that is dCij(t)/dt= 0, where spatial coordinates are 
denoted by (i, j). We assume that the input is nor- 
malized according to e < C < 1, with e is chosen 
such that < e < minj i j{£jj}. Let P denote the 
membrane potential of the photoreceptor, which is 
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FIG. 2: Artifacts with a luminance ramp. Left im- 
age: The input C% , a luminance step with a superim- 
posed luminance ramp (increasing linearly from left to 
the right). Middle image: With the default value 7 = 1.5 
in equation |2l the adaptated image is correctly rendered 
and hardly distinguishable from the input. Right image: 
Setting 7 = causes the appearance of ripple artifacts 
in the adaptated image. All results are shown at t — 250 
iterations. 



FIG. 1: Model sketch. A luminance distribution is 
subjected to a divisive "gain control" stage #1 (S(t), 
equation 13 • At this stage, inhibition of luminance C 
takes place as a function of increasing photoreceptor po- 
tential P. The second gain control stage G(t) can either 
amplify the signal S(t) or attenuate it (equations [HJ 
andHJl. Amplification of S(t) occurs if the membrane 
potential P falls below a threshold O, and attenuation 
for P > Q (equation . Both "gain control" stages 
interact multiplicatively (denoted by the symbol 
equation |5J before providing excitatory input into the 
photoreceptor's membrane potential (symbol "P", equa- 
tion HJ. The photoreceptor potential in turn feeds back 
into both of the "gain control" stages. At the same time, 
the photoreceptor potential represents the output of our 
model. 



assumed to obey the equation (the symbols gi ea k, 
9exc{t), and V exc are defined below) 



dP{t) 
dt 



-gieakP(t) + gexc(t) [V exc - P(t)} (1) 



An instance of the last equation holds for each posi- 
tion hence P = Pij(t) (in what follows indices 
were dropped for brevity). The excitatory satura- 
tion point (or reversal potential) is defined by V exc , 
and the leakage (or passive) conductance is defined 
by gieak (note that V exc represents an asymptote for 
P) . Both of the last constants are equal for all pho- 
toreceptor cells. The default simulation parameters, 
as well as further simulation details, can be found in 
table[n] Notice that photoreceptors in fact hyperpo- 
larize in response to light (c.f. section [nj, whereas 
the last equation makes a contrary assumption. This 
assumption, however, implies no loss of generality, 
since the model can equivalently be re-formulated 
such that it hyperpolarizes with increasing intensity 
levels. 




= 1 (oorrslsnl) f - [no dwrswe gain) fl^i) - fl (conslan1) 



FIG. 3: Artifacts. The results shown in this figure 
should be compared with figure Q Left image: Setting 
the amplification constant to G(t) = 1 in equation [5] di- 
minishes adaptation (i.e., low luminance values are not 
pushed that high). Notice that in this case dynamical 
switching is made inoperative. Middle image: Setting 
7 = in equation has no effect on the natural images 
we have tested, but causes strong ripple artifacts with lu- 
minance ramps as demonstrated in figure|2] Right image: 
Using a constant threshold Q(t) — Oo = 0.25 in equa- 
tion leads to strong saturation (or over-adaptation). 
All results are shown at t = 250 iterations. 



Excitatory input to the photoreceptor potential is 
given by the conductance g exc = g exc ,ijif)i which is 
defined by 



5 (t) = G(t) ■ S(t) 



(2) 



where the process G = Gy(t) interacts multiplica- 
tively with the light-induced signal S = S$ (t) (such 
interaction was previously referred to as mass action 
or gating mechanism, see HJ). For the signal S, we 
assume that its efficiency for driving the photore- 
ceptor's potential diminishes with increasing poten- 
tial P: 



S(t) = 



1 + 7- P(t) 



(3) 



4 



The last equation in fact establishes a feedback 
mechanism which allows the photoreceptor to reg- 
ulate the strength of its own excitatory input. In 
addition, the excitatory drive of the photoreceptor 
is also a decreasing function of increasing potential 
P(t) by virtue of the term "(V exc — P)" (the driv- 
ing potential) in equation ^ Notice that if g exc 
was constant and sufficiently high, the driving po- 
tential would make P(t) saturate at V exc (i.e., V exc 
is asymptotically approached). Therefore, both the 
excitatory input g exc , and the driving potential, de- 
crease as P(t) grows. The motivation for including 
equation in our model was to eliminate ripple ar- 
tifacts seen with luminance ramps (figure EJ). With 
"normal" natural images, those artifacts did not ap- 
pear to be a major nuisance (figure^ see also section 

Eg. 

The process G(t) implements an amplification mech- 
anism as follows: 



n- 



d G(t) 
dt 



-G + S(t-t ) 



(4) 



where 5(t — to) is the Dirac function which serves 
to impose the initial condition G(t = to) = 1, given 
that G{t) = V t < t (notice that G = G l3 {t) as 
usual). Simulations are assumed to start at t = 0. 
By virtue of the index k e {1,2} associated with 
the time constant , the last equation describes two 
distinct processes. These processes are characterized 
by t\ > (making G decay with time), and T2 < 
(leading to an increase of G with time). The last 
equation thus implements what we dubbed a "dy- 
namically switching gain control" . But who or what 
is switching G on (i.e. making it increase with |t"2 |) 
or off (i.e., making it decrease with n)? The one or 
the other process is invoked depending on whether 
P exceeds a threshold 8 or not: 



k = 1 if P{t) > 0(f) 
k = 2 otherwise 



(5) 



This means that if the outer segment potential P is 
below the threshold 0, its input g exc (t) is amplified 
via equation [3] The amplification mechanism acts 
to diminish the integration time of luminance sig- 
nals until reaching the threshold 0, especially low- 
intensity signals. Once the threshold is exceeded, 
amplification is switched off (figure|SJ. In fact, G de- 
cays rapidly then in order to avoid driving the outer 
segment potential into saturation (which neverthe- 
less may occur at sufficiently high intensity values). 
With ineffective dynamical switching G = const, 
adaptation is severely deteriorated (figure|3 first im- 
age). Mathematically, the dynamic switching mech- 
anism avoids an unbounded growth of G. 
Amplification proceeds until P crosses a threshold. 



The threshold, however, is not fixed, but is rather 
represented by a slowly decaying process on its own: 

re^p- = -e(t) + Q -5(t-t ) (6) 

The Dirac function 8(t — t ) establishes the initial 
condition 0(t = t ) = O , if ®(t) = V t < t (no- 
tice that = 0g(t)). We like to emphasize that 
the threshold is not supposed to represent a fir- 
ing threshold for the photoreceptor. It rather serves 
to implement the dynamic switching behavior for 
turning the signal amplification on or off. The mo- 
tivation for including a dynamical threshold in our 
model was the elimination of artifactual contrasts in- 
version effects, and will be explained in more detail 
in section HVI Furthermore, if a constant threshold 
were chosen, over- adaptation would occur (figure |21 
last image). 

Our simulations were evaluated at the moment when 
Pij > Qij V This is, however, not a steady- 

state, because the outer segment potential continues 
to decay with gi ea k ■ The results which are presented 
in the figures 151 to 1 1 01 therefore show snapshots of the 
outer segment potential at exactly the moment when 
the last potential value Pij(t) exceeded the thresh- 
old 0(t) (i.e., corresponds to the position with 
the lowest intensity value in the input). One may 
ask why we gave preference to a dynamical formu- 
lation of our model over steady-state equations. In- 
tuitively, steady-state solutions cannot capture the 
full behavior revealed by the model. For example, 
the steady-state solution (as defined by dQ/dt = 0) 
of the last equation is zero, and, depending on k, the 
steady-state solution of equation0]is infinity (k = 2) 
or zero (A; = 1). 



IV. DESCRIPTION OF THE ADAPTATION 
DYNAMICS 

What does the adaptation dynamics defined by 
equations to [fj] look like? The process obviously 
integrates the activity generated by an input im- 
age C, via the photoreceptor membrane potential P. 
The integration proceeds until P exceeds the thresh- 
old 0. At this point the integration process decel- 
erates exponentially with a time constant t± > 0, 
since the corresponding solution to equation 0] is 
0(<) = exp(— t/ri). The dynamics of P is shown 
in figure^! luminance values that vary over 5 orders 
of magnitude are mapped onto roughly two orders 
of output magnitude in a way that contrast rela- 
tionships of the input are preserved. Moreover, the 
process converges rather fast. Even for the smallest 
input intensities, convergence is reached at about 
200 iterations. This fastness is a consequence of the 
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parameter 


value 


equation 


description 


Qleak 


0.05 


m 


leakage conductance 


Vexc 


1 


m 


synaptic battery 


i 


1.5 


m 


divisive gain 


Tl 


0.7213 


□ 


damping time constant 


T'2 


-40.4979 


m 


amplification time constant 


&o 


0.25 


m 


initial threshold value 


re 


39.4949 


□ 


threshold decay time constant 



TABLE II: Simulation details. The table is 
self-explanatory. For the integration of equa- 
tion a fourth-order Runge-Kutta scheme was 
used with an integration time step of 0.01. 
The remaining differential equations were inte- 
grated with Euler's method with an integration 
time step of one. Notice that the integration 
step sizes were not adjusted to match physio- 
logical time scales. 
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FIG. 4: Photoreceptor potential. The photorecep- 
tor potential P (equation is plotted as a function 
of time (t = to 250 iterations) and input intensity 
(C £ {10 ,10 _1 ,...,10" 5 }. The photoreceptor ampli- 
tude is color-coded (colorbar). "Convergence" occurs 
when the photoreceptor potential P exceeds a threshold 
<d, and corresponds to the area over the diagonal line. 
The minimum integration time is delinated by the hor- 
izontal line at the bottom. With decreasing luminance, 
one observes an increase in integration time until "con- 
vergence" is reached (as illustrated by the red arrows 
pointing to the plateau). A similar increase in integra- 
tion time with decreasing stimulus intensity levels is also 
known from the retina, and is expressed as Bloch's law of 
temporal integration. Bloch's law relates the threshold 
for seeing a stimulus to stimulus duration (i.e. integra- 
tion time) and stimulus intensity: the product of stim- 
ulus duration and stimulus intensity equals a constant 
within a so-called critical time window. Bloch's law is 
especially prominent for scotopic vision. 



dynamic switching process, which increases signal 
amplification G until P exceeds 9 (doing so reduces 
the integration time especially for weak luminance 
signals). Since this process (equation^) per se would 
grow in an unbounded fashion, one may question its 
physiological plausibility. But as long as e > 0, or 
dynamically varying noise is present in the model, 



FIG. 5: Dynamics of the "switching" gain con- 
trol. The same as in figure [I] but here the dynamics of 
the signal amplification variable G(t,C) (equation is 
visualized. The bright (dark) area on the bottom (top) 
indicates where the gain control is switched on (off) . No- 
tice that the switching occurs rather fast around the red 
area. The switching area resembles a blurred line - com- 
pare it to the diagonal line delineating the convergence 
plateau in figure ^] 



eventually all luminance values reach threshold in 
finite time, and as a consequence G (equation |3J 
switches from amplification to attenuation. This is 
to say that for k — 2 the process G is bounded 
mathematically from above. Furthermore, numer- 
ical experiments demonstrate that G does not adopt 
excessively high values (see figure : 42j. 
Nevertheless, a suitably parameterized and asymp- 
totically bounded process for substituting G, rather 
than a sharply cut exponential (as it is implemented 
by the equations 0] and EJ), would perhaps better 
reflect physiological reality - but for the moment we 
set aside plausible functions to keep the model con- 
cise. 

Why should the threshold drop with time? Imag- 
ine that we fix to some constant value. In that 
case, all luminance values are integrated until they 
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the ripples are absent from the input, c.f. figure|3J). 
Those artifacts are counteracted by the additional 
gain control mechanism (equation [3J . Its net effect 
is to continuously decrease the integration step size 
for equation ^ as the potential P grows. This effect 
gets especially prominent for high luminance values 
(see figure El- 
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FIG. 6: Input signal. The same as in figure but 
here the dynamics of the input signal S(t,C) is visualized 
(equation EJ. 



all reach the same threshold. This means that the 
integration process would establish a common level 
for bright and dark luminance values, what in the 
best of all cases would lead to a strong reduction of 
contrasts with respect to the input (hgure|3 last im- 
age). But there is yet another, more technical point, 
to this. 

Consider a pair of luminance values, one brighter 
than the other. Since the integration process pro- 
ceeds with fixed time steps (and exponentially in- 
creasing gain), we may choose both luminance val- 
ues such that they exceed the fixed threshold in a 
way that the previously dark luminance value leaves 
more super-threshold activity than the bright value 
(the brighter value must have exceeded threshold 
at some former time step, and thus its activity P 
already has decayed somewhat due to the passive 
leakage conductance gi ea k hi equation ^) . In other 
words, when decoding the photoreceptor potential 
P, the dark value would suddenly appear brighter 
than the original bright value. Such "contrast in- 
version" artifacts are avoided with a threshold that 
decreases with time. Thus, the dynamic threshold 
process (equation^) acts to preserve contrast polar- 
ities (notice that the threshold process asymptoti- 
cally approaches zero). 

Yet another type of artifact may emerge as a con- 
sequence of the exponentially increasing amplifica- 
tion signal G, most likely due to amplification of nu- 
merical noise while integrating the differential equa- 
tions. With certain luminance distributions, espe- 
cially with luminance ramps, step- like or ripple- like 
structures may appear when P is read out (of course 



FIG. 7: Tiled Lena image. The original Lena image 
(with luminance values between and 1, see figure nTt 
was subdivided into four tiles, and tiles were multiplied 
with 10°, KT 1 , 1CT 2 , and 1CT 3 , respectively. In the in- 
put (leftmost image) , both of the lower tiles are displayed 
in black. The order of magnitude of the corresponding 
luminance range is indicated with the black tiles. 
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FIG. 8: MIT image. The first image shows the input 
image, with luminance values originally varying from 
to 255. The input image was normalized such that the 
maximum intensity value was 1, and the minimum 0. 
Subsequently, all zero luminance values were substituted 
by e — (l/255)/2. The second image shows the result 
obtained with the method described in [HIU (500 iter- 
ations) . The last image was obtained with our approach 
(150 iterations; convergence occurred within simulation 
time). Both results show the cone's membrane potential. 



V. RESULTS OF NUMERICAL 
EXPERIMENTS 

What should one expect from a "good" adaptation 
mechanism? It should map luminance values, which 
can be distributed over several orders of magnitude, 
onto a fixed target range of, say, one or two orders 
of magnitude. In this way, images with a high dy- 
namic range could be visualized with a normal com- 
puter monitor. If we tried a direct visualization of 
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a high dynamic range image without applying any 
adaptation, we could just see the luminance patterns 
of the first one or two orders of magnitude, while all 
smaller luminance values would be displayed in black 
(see figure notice that the optic nerve has a simi- 
lar transmission bandwidth) . Additionally, a "good" 
adaptation mechanism should leave unchanged an 
input image which does only vary over one or two 
orders of magnitude. Or at least leave such an im- 
age as far unchanged as possible. Contrast strength 
should ideally be preserved. Put another way, com- 
pression effects that are introduced by the adapta- 
tion mechanism should be minimized. 
We compare the results of our mechanism with 
one proposed by 0, |2£j (subsequently denoted by 

"g&h") m. 

In order to assure that, at some time, P(t)y > O(t) 
at all positions (i, j), zero values of the original lu- 
minance distribution were substituted by the half 
of the second smallest luminance value, that is e = 
0.5 *min{Cij : > 0}, if not otherwise stated. We 
used standard benchmark images of size 256 x 256 
pixels as inputs C 

Figure [5] show the results with the MIT image, 
where the result obtained with our method is slightly 
less saturated than the one obtained with G&H's 
method. 

In order to better explore the performance of the two 
methods, we superimposed the original test images 
with artificially generated illumination patterns. In 
figure El the MIT- image was multiplied with a lumi- 
nance ramp to simulate an illumination gradient. In 
the latter case, the result from G&H is less saturated 
than ours. 

In figure the original image (shown in figure [TT\ 
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FIG. 9: MIT image with overlying luminance 
ramp. The original MIT image (see figure |SJ was mul- 
tiplied with a luminance ramp which linearly increases 
from left (intensity 0) to the right (intensity 1). 

was subdivided in four "tiles" , where within each tile 
luminance values vary over a different order of mag- 
nitude. This test image mimics a situation where 
the range of luminance values within a scene varies 
over four orders of magnitude. Both methods push 
luminance values sufficiently high such that details 
in the darkest tile are rendered visible (where our 



method yields an overall more brighter result - and 
hence the darkest patch is better visible). Thus, four 
orders of magnitude of input range are mapped onto 
two orders of magnitude available for visualization, a 
situation that is similar to situations which are met 
by the retina. 

In the last example, we created an artificial high- 
dynamic range image (figure I1U[) from the origi- 
nal "Peppers" image ( figure ITTJl . In this case, our 
method produces a slightly brighter result compared 
with G&H: the result generated with G&H's method 
has harder contrasts. 
We conducted further simulations where we set 
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FIG. 10: Power-law-stretched Peppers image. Lu- 
minance values of the original Peppers image (see fig- 
ure 111! were raised to the power of 4 to create a high 
dynamic range image. 




FIG. 11: Original "Lena" and "Peppers" image. 

These images are shown for comparing them with the 
results presented in figure 171 and 1101 respectively. 

Cy <— Py after convergence, and re-started the sim- 
ulation. The results did not change, indicating that 
the model's state after converging the first time al- 
ready corresponds to a steady-state solution. 

VI. MODEL BEHAVIOR WITH 
PARAMETER CHANGES 

The parameters of our model can be tuned 
according to the expected numerical range of 
luminance values. In this way, compression effects 
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in the output are reduced, what can lead to the 
generation of visually more pleasing results. 
Increasing the value of 7 (table |nj equation |2l 
reduces the overall compression of the input at the 
cost of low-intensity regions. This is to say that 
low-intensity regions will appear darker, and regions 
with higher intensities will be rendered with some- 
what improved contrasts. A similar effect results, 
albeit more intense, when increasing the threshold 
decay time constant tq (equation HaJ) . Decreasing 
the initial threshold value Oo (equation EJ) will 
slightly increase overall brightness and compression, 
respectively. The model behavior is quite robust 
against changes in the damping time constant t\, 
since this mechanism is backed up by the signal gain 
control stage (equation . Nevertheless, variations 
in the value of the amplification time constant T2 
bears strongly on the results: a decrease improves 
greatly the adaptation behavior, but if ti is set too 
low artifacts may occur, such as contrast polarities 
being reversed with respect to the input. On the 
other hand, if T2 — > 00, no adaptation at all takes 
place. In future versions of our approach this 
influential parameter could be set automatically as 
a spatially varying function of the structures in the 
input image. 



VII. "THIS THING CALLED EPSILON"... 

As it turned out, a "smart" choice of e can even 
improve the contrasts in the visualization of the re- 
sults. Because for displaying, each image is normal- 
ized to occupy the full range of available gray levels, 
if e is too small with respect to the second small- 
est luminance value, it gets not sufficiently pushed 
by the adaptation process, such that in the adapted 
image the difference between the smallest and the 
second smallest value is too big. As a consequence, 
many of the darker gray levels are not used (if we 
assume a linear mapping of activity to gray lev- 
els), what leaves less gray levels for displaying the 
other (higher) luminance values. Hence, the con- 
trasts in the displayed image will be reduced. Ide- 
ally, e should depend in some way on how dark the 
input image is perceived by a human observer. Find- 
ing an adequate function that automatically sets the 
value of e would be an interesting topic for future re- 
search. 



VIII. DISCUSSION AND CONCLUSIONS 

We presented a novel theory about the adapta- 
tional mechanisms in retinal photoreceptors. Our 



theory is abstract in the sense that we did not at- 
tempt to identify model stages with components of 
the phototransduction cascade (as outlined in sec- 
tion^. Nevertheless, one is tempted to draw corre- 
sponding parallels between our model and physiolog- 
ical data. In the transduction cascade there are (at 
least) two sites of amplification: the serial activation 
of transducins by the active form of rhodopsin Rh* , 
and the hydrolysis of cGMP by phosphodiesterase. 
An amplification of the signal takes also place in 
our model by virtue of G in equation 01 Further- 
more, Ca 2+ constitutes a messenger for adaptation. 
In contrast, there is no corresponding variable for 
describing the concentration of Ca 2+ in our model. 
Nevertheless, the membrane potential P subserves 
two different purposes. First, it corresponds to the 
output of the photoreceptor. Second, it constitutes a 
feedback signal that acts to control signal amplifica- 
tion - and hence the adaptation process. As Ca 2+ is 
known to be linearly related to the membrane poten- 
tial, it seems reasonable to consider P as a lumped- 
together description for both the membrane poten- 
tial and the Ca 2+ concentration. 
Indeed, one can draw further parallels. In our model, 
signal amplification stops as soon as the membrane 
potential exceeds a threshold, in order to counteract 
saturation effects (equation[SJ). This process is remi- 
niscent on the binding of arrestin to phosphorylated 
Rh*, leading to a complete inactivation of the pho- 
topigment, and thus to a ceasing of the transduction 
cascade. 

In our model, there is yet another way to counteract 
saturation effects, by means of the divisive inhibition 
stage (equation[3J). This process can be compared to 
the interaction of Ca 2+ with the visual cycle, which 
causes an acceleration of the rate of Rh* phospho- 
rylation [36l l37l l3S| . This interaction is brought 
about by the Ca 2+ -binding protein recoverin, and 
decreases the lifetime of Rh*. As a consequence, less 
cGMP will be hydrolyzed upon absorption of a pho- 
ton HJ. 

On the technical side, computer simulations demon- 
strated that our approach is on a par with a rece ntly 
proposed model from Grossberg and Hong [l^. |20| 
("G&H"). However, several crucial differences exist 
between their approach and ours. 
First and above all, the critical stage for adapta- 
tion in G&H's approach consists of the feedback pro- 
vided by electrically coupled horizontal cells. Light 
adaptation through the outer segment can be de- 
coupled from the actual adaptation dynamics, and 
hence may be considered as a pre-processing step in 
their model. 

Remarkably, our approach achieves similar adapta- 
tion results without incorporating the horizontal- 
to-cone feedback loop. This prediction is consis- 
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tent with physiological data, as cone photorecep- 
tors can decrease their sensitivity over about 8 log 
units of background intensity [lj|- Moreover, feed- 
back from horizontal cells may even further im- 
prove adaptation. Since we have seen, on the other 
hand, that contrasts are reduced as a consequence 
of the dynamic range compression, one may specu- 
late that feedback from horizontal may also com- 
pensate for this effect, by re-enhancing contrasts. 
Notice that contrast enhancement is tantamount to 
center-surround interactions. Because adjacent hor- 
izontal cells of the same type are fused by gap junc- 
tions, their feedback will influence the membrane 
potential of neighboring cones within some radius 
of the actually stimulated photoreceptor. In this 
way, the antagonistic receptive field structure is cre- 
ated in bipolar cells. But then bipolar cells repre- 
sent a contrast-enhanced signal of the photorecep- 
tors. Therefore, neurophysiological data are consis- 
tent with our ideas. 

Both models are similar complex with respect to pa- 
rameter spaces. G&H's approach has some 10 pa- 
rameters, whereas ours has 7 (plus the e). Although 
we did not carry out a detailed analysis of com- 
putational complexity, the respective model struc- 
tures suggest that the G&H model is computation- 
ally more demanding. The latter fact seemed to 
be confirmed with our simulations on a serial com- 
puter, where our model converged in a fraction of 
the time that was necessary to achieve comparable 
results with the G&H model 0] 
Similar to the G&H model, another approach [3jJ is 
also motivated by the observation that strong con- 
trasts usually indicate reflectance changes in natu- 
ral scenes, as opposed to intensity variantions due 
to changes in illumination. The approach from |3^| , 



however, has no stage for luminance adaptation, and 
only computes an "anisotropically-like" smoothed 
version of the image, which is used for exerting di- 
visive gain control directly on intensity values (c.f. 
table P). The lateral connectivity between cells that 
form the diffusion layer is controlled by inverse We- 
ber contrasts. Hence, both strong and weak con- 
trasts in the original image may affect the degree of 
smoothing. Simulation results obtained with our im- 
plementation of Gross' and Brajovic's approach re- 
vealed strong boundary enhancement if tuned such 
that the adaptation was comparable to the other two 
methods. This suggests that the signal transduction 
characteristics of Gross' and Brajovic's approach is 
high-pass. 

Our model, perhaps with different parameter values, 
should as well be useful for displaying high-dynamic 
range images, or synthetic aperture radar images. 
This is a topic that will be pursued with future re- 
search. Further interesting questions address the in- 
corporation of feedback from horizontal cells, and 
possibly of reset mechanisms for the threshold pro- 
cess, in order to extend our model's processing ca- 
pacities to image sequences. 
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